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In the qualitative theory of ordinary differential equations, we can find many papers 
whose objective is the classification of all the possible topological phase portraits of a given 

co ■ family ° f differential system ' Most . ° f the studies rely ° n systems with real parameters and 

^| ■ the study consists of outlining their phase portraits by finding out some conditions on the 

parameters. Here, we studied a susceptible-infected-susceptible (SIS) model described by the 
differential system x = —bxy — mx + cy + mk, y = bxy — (m + c)y, where 6, c, k, in are real 
parameters with b ^ 0. m ^ [5]. Such system describes an infectious disease from which 

■ infected people recover with immunity against reinfection. The integr ability of such system 
has already been studied by Nucci and Leach [S] and Llibre and Vails [5] . We found out two 
different topological classes of phase portraits. 

O ' 

■ Key-words: SIS epidemic model, global phase portrait, endemic and disease-free steady 
' states. 



1 Introduction 

There is a long time that scientists have been curious about the interaction between 
portions of a population with particular characteristics, e.g. prey and predator interaction 
[2]. In 1838, Verhulst [TJ] proposed the study of the population growth by means of the so 
called logistic equation (see also [2]). In contrast, this equation has had other applications, 
for instance in the study of the spread and the evolution of diseases. The means of how 
diseases spread has stimulated the interest mainly within researchers of the biological field 

ma. 

As suggested above, the application of mathematical tools in different areas of sciences has 
been frequent and has had a great impact on the scientific and/or experimental conclusions. 
For example, Lotka (1925) and Volterra (1926) described independently the dynamics present 
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in biological interactions between two species using ordinary differential equations (ODE). 
Their model is now known as the Lotka-Volterra equations and they are the most studied 
equations in the qualitative theory of ODE (see (2J[3] for more details). 

Another approach of the ODE is their application in the study of infectious diseases. 
According to Levin [5], this study represents one of the oldest and richest areas in mathe- 
matical biology. Besides, continuous models composed by ODE have formed a large part of 
the traditional mathematical epidemiology literature, mainly because mathematicians have 
been attracted by applying the ODE's tools, such as their qualitative theory, to the study of 
infectious diseases in the attempt of using mathematics to contribute positively to the science 
field and because the mathematical models become indispensable to inform decision-making. 

Particularly, we consider the system of first-order ODE 

x = —bxy — mx + cy + mk, 
y = bxy - (m + c)y, 

where x and y represent, respectively, the portion of the population that has been susceptible 
to the infection and those who have already been infected. System ([T]) is a particular case 
of the class of classical systems known as susceptible-infected-susceptible (SIS) models, in- 
troduced by Kermack and McKendrick [3] and studied by Brauer [3J , who has assumed that 
recovery from the nonfatal infective disease does not yield immunity. In system ((T]), k is the 
population size (susceptible people plus infected ones), mk is the constant number of births, 
m is the proportional death rate, b is the infectivity coefficient of the typical Lotka-Volterra 
interaction term and c is the recovery coefficient. As system ((TJ) is assumed to be nonfatal, 
the standard term removing dead infected people —ay in [3] is omitted. As usual in the lit- 
erature, all the critical points of system (JlJ will henceforth be called (endemic) steady states 
(e.g. see [13]). In the rest of the paper we will study the phase portraits of the differential 
system (TTJ) with bm ^ 0. Note that if b — 0, then system ([T]) becomes linear, and if m = 0, 
then system ([T]) satisfies that x + y — 0. These last two cases are trivial and non interesting 
from a biological point of view. 

Much has been studied on SIS models. The great part of the studies present only local 
stability results for which many mathematical tools are available. In contrast, global studies 
of these models are very limited due to the lack of applicable theories. The main mathe- 
matical tool which has been used for this purpose is the Lyapunov function. The principal 
limitation is that such functions are not the same for all types of systems (see [TTJ [T^] and 
references there listed for further details). 

Besides the stability of SIS models, their integrability has also been studied. For example, 
Nucci and Leach [5] have demonstrated that (UJ is integrable using the Painleve test. Later, 
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Llibre and Vails |6J have proved that system ([TJ) is Darboux integrable, and they have shown 
the explicit expression of its first integral and all its invariant algebraic curves. 

Alternatively, the attempt of outlining the global phase portraits of ODE systems is a 
possible way to determine their global behavior. 

Here, the purpose was to classify all the topological classes of the global phase portraits 
of system ([1]) using some information in [6J . 

The main result in this paper is the following: 

Theorem 1.1. The phase portrait on the Poincare disc of system ((T|) is topologically equiva- 
lent to one of the two phase portraits shown in Figure [7J modulo reversibility. 



The paper is organized as follows. Basic definitions and stated results about singular 
points on the plane and their classification are given in Section [SJ The reader who is familiar 
to these concepts can skip this section with no worries of missing new information. In Section 
[3l the local analysis of the finite and infinite singular points of system (JXJ) is done, and so is 
the existence of at least three invariant straight lines for such system. Section 0] collects all 
the information provided in the previous sections and proves Theorem ll.il Finally, Section 
[5] concludes and discusses the results from the biological point of view. 

2 Basic results 

In this section, we introduce some definitions and stated theorems which shall be used in 
the next section for the analysis of the local phase portraits of finite and infinite singularities 
of system (TTJ) . 




(a) 



(b) 



Figure 1: Phase portraits in the Poincare disc of system ([!]) 
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2.1 Singular points 

We denote by X(x, y) = (P(x, y), Q(x, y)) a vector field in R 2 , associated to the system 



x = P(x,y) 
y = Q(x,y) 



(2) 



where P and Q are real polynomials in the variables x and y of degree at most neM. We 
recall that the degree n is defined by the maximum of deg(P) and deg(Q). We denote by 
P n (M. 2 ) the set of all polynomial vector fields on R 2 of the form ©. 

Definition 2.1. A point pel 2 is said to be a singular point, or a singularity, of the 

vector field X = (P, Q), if P(p) = Q(p) = 0. 

We denote by P x the partial derivative of P with respect to x. Suppose that p G R 2 is 
a singular point of X and that P and Q are analytic functions in a neighborhood of p. Let 
S = P x {p)Q y {p) - Py(p)Q x (p) and r = P x (p) + Q v {p). 

Definition 2.2. The singular point p is said to be non-degenerate ifS^O. 

As a consequence, p is an isolated singular point. Furthermore, p is a saddle if S < 0, a 
node if r 2 — Ad > (stable if r < 0, unstable if r > 0), a focus if r 2 — A5 < (stable if r < 0, 
unstable if r > 0), and either a lueafc focus or a center if r = and <5 > [7]. 

Definition 2.3. T/ie singular point p is said to be hyperbolic if the two eigenvalues of the 
Jacobian matrix J X(p) have nonzero real part. 

Hence, the hyperbolic singularities are the non-degenerate ones, except the weak foci and 
the centers. 

Definition 2.4. The degenerate singular point p (i.e. 8 = 0), with r ^ 0, is called semi- 
hyperbolic. 

Again p is isolated in the set of all singular points. Specifically, the next theorem summa- 
rizes the results on semi-hyperbolic singularities that shall be used later. For more details, 
see Theorem 2.19 of [7]. 

Theorem 2.5. Let (0,0) be an isolated singular point of the vector field X given by 



where A and B are analytic in a neighborhood of the origin starting with at least degree 2 
in the variables x and y. Let y — f(x) be the solution of the equation Xy + B(x,y) = in 



x 



y 



A(x,y), 

Xy + B(x,y) 



(3) 
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a neighborhood of the point (0,0), and suppose that the function g(x) = A(x,f(x)) has the 
expression g{x) — ax a + o(x a ), where a > 2 and a ^ 0. So, when a is odd, then (0,0) is 
either an unstable node, or a saddle, depending if a > 0, or a < 0, respectively. In the case of 
the saddle, the separatrices are tangent to the x-axis. If a is even, the (0, 0) is a saddle-node, 
i.e. the singular point is formed by the union of two hyperbolic sectors with one parabolic 
sector. The stable separatrix is tangent to the positive (respectively, negative) x-axis at (0,0) 
according to a < (respectively, a > 0). The two unstable separatrices are tangent to the 
y-axis at (0, 0). 

Definition 2.6. The singular points which are non-degenerate or semi-hyperbolic are called 
elementary. 

Finally, we present two other concepts that shall be used in the next sections. 

Definition 2.7. A function H : U C K 2 — » R of class C 1 is called a first integral of system 
(|2"j) if U is an open and dense subset of R 2 and H is constant over the solutions of system 
(T2J) contained in U, i.e. 

OH BH n 
ox ay 

on the points of U . If such H exists, then we say system Q is integrable. 

Definition 2.8. Let C[x, y] denotes the ring of complex polynomials in the variables x and 
y. The polynomial f(x,y) £ <C[x,y] \ C is invariant under system © if 

P^ + Q d / = Kf, 

ox ay 

for some K £ C[x,y] called the cofactor of f. 

We say the algebraic curve f(x,y) = is invariant, if f{x,y) is invariant under system 

©• 

2.2 Poincare compactification 

Let X £ P„(R 2 ) be a planar polynomial vector field of degree n. The Poincare compacti- 
fied vector field 7r(X) corresponding to X is an analytic vector field induced on § 2 as follows 
(for more details, see [7]). 

Let S 2 = {y = (yi, j/2, 2/3) S R 3 ; y\ +y% + y| = 1} and T y S 2 be the tangent plane to § 2 at 
point y. Identify R 2 with T^ 01 ^S 2 and consider the central projection / : T( ,o,i)§ 2 — ^ § 2 - 
This map defines two copies of X on S 2 , one in the northern hemisphere and the other in 
the southern hemisphere. Denote by X' the vector field D/ o X defined on S 2 except on its 
equador S 1 = {y £ § 2 ; y 3 = 0}. Then, S 1 is identified to the infinity of R 2 . 
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In order to extend X' to a vector field on S 2 , including S 1 , X should satisfy suitable 
conditions. In the case that X G P n (M. 2 ), tt(X) is the only analytic extension of y$~ X' 
to § 2 . On § 2 \ S 1 there exist two symmetric copies of X, and knowing the behavior of 
ir(X) around S 1 , one can conclude the behavior of X in a neighborhood of the infinity. The 
Poincare compactification has the property that S 1 is invariant under the flow of tt(X). The 
projection of the closed northern hemisphere of § 2 on 2/3 = under (yi, y 2 , 2/3) 1— > (yi, 2/2) is 
called the Poincare disc, and it is denoted by D 2 . 

Definition 2.9. Two polynomial vector fields X and Y on M 2 are topologically equivalent 

if there exists a homeomorphism on S 2 preserving the infinity S 1 carrying orbits of the flow 
induced by n(X) into orbits of the flow induced by ir(Y). 

As S 2 is a differentiable manifold, for computing the expression for n(X), consider six 
local charts Ui = {y G S 2 ; yi > 0} and Vi = {y G S 2 ; yi < 0}, where i = 1,2,3, and the 
diffeomorphisms Fi : Ui — >• M 2 and Gi : Vi — > M 2 , for i — 1,2,3, which are the inverses 
of the central projections from the tangent planes at the points (1,0,0), (—1,0,0), (0,1,0), 
(0,-1,0), (0,0,1) and (0,0,-1), respectively. Denote by z = (u,v) the value of Fi(y) and 
Gi(y), for any i = 1, 2, 3 (so z represents different things according to the local charts under 
consideration). Hence, after some easy computations, tt(X) has the following expressions: 

W z )( Q (l^)-rt(l») t - v p(l»)) inU u (4) 

\ \V V J \V V J \V V J J 

v»A(z)(p(^)-uq(^),-vq(?-,1)) mU 2 , (5) 
\ \v v) \v vj \v vJJ 

A(z)(P(u,v),Q(u,v)) in U s , (6) 

where A(z) = (u 2 + v 2 + l) _ ( n_1 )/ 2 . The expression for Vi is the same as that for Ui except 
for a multiplicative factor (— 1)"~ 1 . In these coordinates for i = 1,2, v = always denotes 
the points of S 1 . 

After finding the infinite singular points, it is necessary to classify them. Among the 
hyperbolic singular points at infinity only nodes and saddles can appear. All the semi- 
hyperbolic singular points can appear at infinity. 

If one of these hyperbolic or semi-hyperbolic singularities at infinity is a (topological) 
saddle, then the straight line {v — 0}, representing the equator of § 2 , is necessarily a stable 
or unstable manifold, or a center manifold (see Figure [5]). 

The same property also holds for semi-hyperbolic singularities of saddle-node type. They 
can hence have their hyperbolic sectors split in two different ways depending on the Jacobian 
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matrix of the system in the charts U\ or U^. The Jacobian matrix can be either 













, or 


: 


: 











with A ^ 0. In the first case we say that the saddle-node is of type SNl and in the second 
case of type SN2. The two cases are represented in Figure [3J The sense of the orbits can also 
be the opposite. 




S SNl SN2 

Figure 2: (S) represents a hyperbolic or semi-hyperbolic saddle on the equator of § 2 ; (SNl) and 
(SN2) represent, respectively, saddle-nodes of types SNl and SN2 of vr(X) in the equator of § 2 



3 Analysis of the system 

Here, we provide a mathematical analysis of system (fTJ) . 

Provided that 6^0, system (TTJ) has only two finite singular points: 
p = ((c + m)/b, (— c + bk — m)/b), usually known as endemic steady state, and q — (fc,0), 
usually known as disease-free steady state. In addition, both finite singular points p and q 
are the same if bk = c. 

The names given above to each steady states are not an accident. In q, the number of 
susceptible individuals is equal to the population size fc, whereas the number of infected 
people is null. On the other hand, the number of susceptible people in p is the recovery 
coefficient plus the death rate divided by the infection coefficient, while the infected ones 
are the rest of the population, which leads to the presence of infected people, since bm ^ 0. 
Finally, note that only non-negative values of x and y are interesting here, because they 
represent the number of individuals. 

First, we start with the analysis of the endemic steady state p. Translating the singular 
point p to the origin in system (TTJ), we obtain 

x = —bkx + cx — my — bxy, 
y = (— c + bk — m)x + bxy, 



(7) 
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which is equivalent to ([T]). The Jacobian matrix of 10 is given by 



S(x,y) 



( 



c + bk — m + by 



c — bk — by 



m — bx 



bx 



) 



which implies that 



5 = (5(0, 0) = (bk — c — m)m and r 



t(0,0) = -bk + c. 



If (bk — c — m)m < 0, then p is a saddle point. On the other hand, if (bk — c — m)m > 0, 
then p is a node point, because r 2 — 46 = (c — bk + 2m) 2 > 0. 

In the case that (bk — c — m)m = 0, or equivalently, m — bk — c, then p is degenerate. 
Indeed, it is the case that both finite singular points are the same, i.e. p = q = (fc, 0). Here, 
system ([7]) becomes 



so that p is a semi-hyperbolic point. By a linear change of coordinates, system ([5]) can be put 
on the form of system @ and, applying Theorem 12. 51 we conclude that p is a saddle-node 
point. 

From the biological point of view, when the death rate m is equal to the portion of the 
population which becomes infected (bk) minus the recovery coefficient, the dynamics around 
the steady states p and q changes and they become only one point which attracts (the node 
part) and repels (the saddle part) the orbits in its neighborhood. 

Now, we analyze the disease-free steady state q. Translating the singular point q to the 
origin in system (fTJ) , we obtain 

x = —rax + cy — bky — bxy, 
y = (—c + bk — m)y + bxy, 

which is equivalent to (TTJ). The Jacobian matrix of © is given by 



x 



mx — my — bxy 



(8) 



V = bxy, 



whose Jacobian matrix at (0, 0) is 
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which implies that 

6 = 5(0, 0) = -(bk - c - m)m and t = r(0, 0) = bk - c - 2m. 

If — (bk — c — m)m < 0, then q is a saddle point. In contrast, if 
— (bk — c — m)m > 0, then q is a node point, because r 2 — AS = (c — bk) 2 > 0. 

The case (bk — c — m)m = 0, or equivalently m = bk — c, has already been studied, and 
p = q is a semi-hyperbolic saddle-node point. 

Finally, we have proved the following: 

Proposition 3.1. Consider system ([1]) with bm ^ and its two finite steady states p and 
q. Then: 

1. If either m > and m > bk — c, or m < and m < bk — c, then p is a saddle and q is 
a node; 

2. If either m > and m < bk — c, or m < and m > bk — c, then p is a node and q is a 
saddle; 

3. If m = bk — c, then p = q is a semi-hyperbolic saddle-node. 

Having classified all the finite singular points, we apply the Poincare compactification to 
study the infinite singularities. 

In the local chart U\, where x — 1/v and y = u/v, we have: 

u = u(b + bu — cv — cuv — kmv 2 ), 

(10) 

v = v(bu + mv — cuv — kmv), 

whose infinite singular points are (0, 0) and (—1, 0), which are a saddle-node of type SN1 (by 
Theorem 12. 5|) and a node, respectively. 

In the local chart U2, where x — u/v and y = 1/v, the system 

u = —bu — bu 2 + cv + cuv + kmv 2 , 

(11) 

v = v(—bu + cv + mv) 

has two infinite singular points (0,0) and (—1,0). The latter one is a node and is the same 
as (—1,0) G U%, while the former one is a saddle-node of type SN1. 
We have just proved the following: 

Proposition 3.2. The infinite singular points of system (TI| are the origin of charts U\, V\, 
U2 and V2, which are saddle-node points of type SN1, and (—1,0), belonging to each of the 
charts U\ and U2, which is a node point. 
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Knowing the local behavior around finite and infinite singular points, another useful tool 
to describe the phase portraits of differential systems is the existence of invariant curves. 
The next result shows system ([T]) has at least three invariant straight lines. 

Proposition 3.3. Let bm ^ 0. System ([T]) has at least three invariant straight lines given 
by fi(x,y) — y and f2{x 1 y) = k — x — y, and additionally fz{x,y) — k — x, if c = bk. 

Proof. By Definition 12.81 we can find K\(x, y) = bx — m — c, Ki{x, y) = —m and K^(x, y) — 
—m — by as the cofactors of fi(x,y), f-2{x,y) and fa(tc,y) (if c = bk), respectively. □ 

4 Main result 

From Propositions 13.11 and 13.21 we get all the information about the local behavior of 
finite and infinite singular points, respectively. Using the continuity of solutions and primary 
definitions and results of ODE (e.g. w-limit sets, existence and uniqueness of solutions, the 
Flow Box Theorem etc. [7j) and the existence of invariant straight lines of system (JTJ stated 
by Proposition 13.31 its global phase portraits can be easily drawn. 

Essentially, we have only two cases. The finite steady state q is the intersection of the 
invariant curves fi(x, y) = f^{x, y) — 0, and the other finite steady state p lies on the curve 
h{x,y)=Q. 

According to items (1) and (2) of Proposition 13.11 p (respectively, q) is a saddle (re- 
spectively, a node) the one way and the other a node (respectively, a saddle). The subtle 
difference here is the position of point p. While q remains on the line {y — 0}, p is in the 
lower part of the Poincare disc the one way and the other in the upper part. It is worth 
mentioning that the four infinite singular points continue to be the same points no matter 
what conditions are being considered. The phase portrait of both cases above is topologically 
equivalent to the one which is shown in Figure [IJa). 

Item (3) of Proposition 13.11 assumes the existence of only one finite singular point, p = q. 
Here, when m — bk — c, both p and q become only one degenerate singularity which bifurcates 
into a saddle-node point. Again, no changes are applied to the infinite singular points. The 
phase portrait of this case is topologically equivalent to the one which is shown in Figure 

raw. 

Finally, Theorem 1 1.1 1 has been proved. 

5 Conclusions and discussions 

The last section proves the existence of only two classes of global phase portraits of the 
quadratic system (|TJ) . In the qualitative theory of ODE it is quite important to know the 
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global behavior of solutions of systems and, in general, this is not an easy task. The most 
frequently used tools for this propose are the study of local behavior along with the (local 
and global stability), integrability, and also the global phase portrait, which was employed 
in the present study. 

In the case represented by Figureflja), it is clear that while the steady state q characterizes 
the presence of only susceptible individuals, p indicates the mutual presence of susceptible 
and infected people. Besides, as q is an asymptotically stable node, the disease seems to 
be controlled and the whole population tends to be healthy but susceptible to be infected 
again. As p is an unstable saddle steady state, it suggests that there is no harmony between 
the number of susceptible people and infected ones, although some of the solutions tend to 
q. indicating the control of the disease. 

In case of Figure QJb), all the solutions tend to q (regarding that x,y > 0), i.e. if 
m = bk — c, the disease is supposed to be controlled and the whole population is inclined to 
be healthy but susceptible to the reinfection. 
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